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ABSTRACT 

Mini-Neptunes and volatile-poor super-Earths coexist on adjacent orbits in proximity to host stars such as Kepler-36 and 
Kepler-11. Several post-formation processes have been proposed for explaining the origin of the compositional diversity between 
neighboring planets: the mass loss via stellar XUV irradiation, degassing of accreted material, and in-situ accumulation of the 
disk gas. Close-in planets are also likely to experience giant impacts during the advanced stage of planet formation. This study 
examines the possibility of transforming volatile-rich super-Earths / mini-Neptunes into volatile-depleted super-Earths through 
giant impacts. We present the results of three-dimensional hydrodynamic simulations of giant impacts in the accretionary and 
disruptive regimes. Target planets are modeled with a three-layered structure composed of an iron core, silicate mantle and 
hydrogen/helium envelope. In the disruptive case, the giant impact can remove most of the H/He atmosphere immediately and 
homogenize the refractory material in the planetary interior. In the accretionary case, the planet is able to retain more than 
half of the original gaseous envelope, while a compositional gradient suppresses efficient heat transfer as the planetary interior 
undergoes double-diffusive convection. After the giant impact, a hot and inflated planet cools and contracts slowly. The extended 
atmosphere enhances the mass loss via both a Parker wind induced by thermal pressure and hydrodynamic escape driven by the 
stellar XUV irradiation. As a result, the entire gaseous envelope is expected to be lost due to the combination of those processes 
in both cases. Based on our results, we propose that Kepler-36b may have been significantly devolatilized by giant impacts, while 
a substantial fraction of Kepler-36c’s atmosphere may remain intact. Furthermore, the stochastic nature of giant impacts may 
account for the observed large dispersion in the mass-radius relationship of close-in super-Earths and mini-Neptunes (at least to 
some extent). 

Subject headings: equation of state — hydrodynamics — planets and satellites: formation — planets and 
satellites: interiors — stars: individual (Kepler-36, Kepler-11) 


1. INTRODUCTION 

Data obtained from the Kepler mission (Borucki et al. 2010) 
indicate an abundant population of close-in planets with radii 
of Rp = 1 - 4/^0 (Petigura et al. 2013; Burke et al. 2014; 
Dressing & Charbonneau 2015). Although the masses of only 
a handful of those planets have been measured with radial- 
velocity or transit-timing variations (TTV) observations, a 
large range of their mean densities Pp indicates that they have 
diverse bulk compositions (e.g. see Fig.4 in Dressing et al. 
2015). These kinematic and structural properties impose con¬ 
straints on theories of planetary origin. 

It has been suggested that these planets may have formed 
through the coagulation of grains and planetesimals at their 
present-day location (Hansen & Murray 2012; Chiang & 
Laughlin 2013). Although all volatile grains in their natal 
disks are sublimated in the proximity of their central stars, 
embryos may have formed prolifically from refractory grains 
(Li et al. submitted) and acquired or retained a modest amount 
of hydrogen/helium atmosphere in situ with relatively low ap¬ 
parent mean densities (Ikoma & Hori 2012). 

Close-in planets are exposed to intense XUV (extreme UV 
and X-ray) irradiation from their host stars. Photoevapora¬ 
tion can significantly modify the structure of their atmosphere 
(e.g., Owen & Wu 2013). One possible cause for the appar¬ 
ent dispersion in pp is the retention efficiency of H/He atmo¬ 
spheres among planets around host stars (Lissauer et al. 2013) 
with different XUV luminosities during their pre-main se¬ 
quence stage. But around some host stars, high-density super- 


Earths and low-density hot-Neptunes with similar Mp and P 
(where Mp is the planetary mass and P is the orbital period) 
coexist. Their pp diversity cannot be simply attributed to the 
XUV flux from their common host stars. 

Around a particular sub-giant star, two planets, Kepler-36b 
and c, were found (Carter et al. 2012) with adjacent orbits (P = 
13.8 and 16.2 days), comparable masses (Mp = 4.46 ±0.3 
and 8.10 ± 0.53 M 0 ), and radii (Rp = 1.48 ± 0.03 and 
3.67 ±0.05 R 0 ). The inferred pp of the inner super-Earth 
(Kepler-36b) matches that of an Earth-like composition and 
that of the outer mini-Neptune (Kepler-36c) is consistent with 
an internal composition of a thick hydrogen/helium (H/He) 
envelope (8.6 ± 1.3 wt. %, assuming an Earth-like core) atop 
a rocky/iron core (Lopez & Fortney 2013). Although a frac¬ 
tion of Kepler 36c’s atmosphere may be due to the degassing 
of accreted planetesimals (Elkins-Tanton & Seager 2008) or 
in-situ the accumulation of a residual disk gas onto its solid 
core (Ikoma & Hori 2012). However, such a large density 
contrast between these two closely-spaced planets is incom¬ 
patible with either of these two accretion processes. 

It is possible that these two planets may once have had a 
similar amount of atmospheric H/He and their compositional 
dichotomy could be caused by their different core masses 
(Lopez & Fortney 2013). Provided that the lost mass is a small 
fraction of Mp, their near 7:6 mean motion resonance (MMR) 
orbits would not be significantly affected by it (Teyssandier 
et al. 2015). However, Kepler 36c appears to have retained its 
atmosphere and Kepler 36b appears to have not, albeit the for- 
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mer has a smaller surface escape velocity and is expected to 
be more vulnerable to photoevaporation than the latter. An¬ 
other system with a similar dichotomy is Kepler 11b and c 
(Lissauer et al. 201 1). 

Another scenario is that protoplanetary embryos formed 
throughout the disk, including the volatile-rich cold outer 
regions, and converged through type I migration to their 
present-day locations in the proximity of their host stars (Pa- 
paloizou & Terquem 2010; Ida et al. 2013). According to 
the conventional accretion scenario, protoplanetary embryos 
emerge through oligarchic growth and experience cohesive 
collisions with embryos with comparable masses (Kokubo & 
Ida 1998). The energy released in such catastrophic events 
intensely heats both the mantle and atmosphere and may in¬ 
duce substantial losses of planetary atmosphere (Schlichting 
2014). The giant impact hypothesis is consistent with the 
chaotic dynamics exhibited in the Kepler-36 system (Deck 
et al. 2012). Hydrodynamic and A-body simulations have re¬ 
vealed that stochastic/convergent migration can lead to close 
encounters between the two planets and physical collisions 
with other embryos (Paardekooper et al. 2013; Quillen et al. 
2013). 

Due to the stochastic nature of these giant impact events, 
structural diversity among members of closely packed, multi¬ 
ple short-period planet systems is the hallmark of oligarchic 
growth, extensive migration, and giant impacts. This diversity 
takes the form of different levels of shock processing, dif¬ 
ferent final material inventories and planetary densities (As- 
phaug 2010). 

In this paper we explore the possibility that structural di¬ 
versity is established by giant impacts during the formation of 
compact multiple-planet systems. Collisions and disruptions 
of differentiated super-Earths with terrestrial compositions 
(Marcus et al. 2009) or water : rock ratios of 1 : 1 (Marcus 
et al. 2010) have been previously simulated with the smoothed 
particle hydrodynamics (SPH) scheme (Reufer et al. 2012). 
However, there have been no systematic simulations of giant 
impacts between super-Earths with H/He atmospheres. 

In §2, we outline the numerical method and the initial and 
boundary conditions used for these simulations. In our mod¬ 
els, we adopt physical parameters similar to those of the Ke¬ 
pler 36 system and scrutinize the collisional origin of the 
compositional diversity between Kepler-36b and c. In §3, 
we present the results of two sets of numerical simulations 
to show that energetic giant impacts can cause strong shocks 
and lead to the ejection of a large amount of the H/He en¬ 
velope. Shock dissipation also heats up the planetary at¬ 
mosphere and interior, elevating the mass-loss rate induced 
by hydrodynamic escape and radiative evaporation (Hebrard 
et al. 2004). We summarize our results and discuss their im¬ 
plications in §4. 

2. COMPUTATIONAL METHOD 

Our three-dimensional hydrodynamic simulations are based 
on the framework of ELASH (Eryxell et al. 2000), an Eu- 
lerian code with an adaptive-mesh refinement (AMR) capa¬ 
bility. This finite-volume scheme has advantages over most 
commonly-used SPH codes in terms of capturing shock waves 
as well as characterizing low-density and sparse regions. Be¬ 
sides, advection of different fiuid species in a grid cell is al¬ 
lowed to obey their own advection equations, which is ideally 
suited for mixing problems. The main purpose of this work 
is to investigate the general aftermath of giant impacts, so our 
study here is on head-on collisions. 


In order to minimize the artifact of artificial diffusivity 
(Tasker et al. 2008), collision simulations are performed in 
a pair of planets’ center-of-mass frame rather than in the iner¬ 
tial reference frame in which they move through the numerical 
grids with a Keplerian speed. The width of the computational 
domain is 1 x 10^^ cm on each side, and we choose an open 
boundary condition for our simulations. A multiple expan¬ 
sion with a large angular number (L^ax = 60) about the center 
of mass of the embryo pair is adopted in order to compute the 
self-gravity with a sufficient angular resolution. 

Most hydrodynamic simulations of giant impacts ignore the 
presence of a central star and its tidal force. In the present 

1 /Q 

context, the Hill radius of a planet, ri{ = a {q/3) (where 
a is the semi-major axis of the planet and q is the mass ratio 
of the planet to its parent star) is only an order of magnitude 
larger than its physical size. Therefore, the tidal force of the 
central star may have a significant effect on the planet’s post¬ 
impact evolution. 

The orbit of the embryo pair around the host star is deter¬ 
mined by solving a two-body problem that treats a central star 
and the embryo pair as two separate point masses (Guillochon 
et al. 2009). Here we consider giant impacts between a target 
planet Mj and an Earth-mass impactor Mi at 0.1 AU separa¬ 
tion from a 1 Mq star. Initially, the target and impactor are 
held apart from their center of mass by a few target’s radii un¬ 
til they have relaxed into a hydrostatic equilibrium subjected 
to each other’s tidal perturbation. They are released with their 
relative velocity gradually increased from rest to the values 
listed in Table 1 . Eurther acceleration is determined by their 
mutual gravity. 

Two simulations of head-on collisions in the accretionary 
and disruptive regimes are performed. In the accretionary 
model 1, Mt = 4.3 i.e. close to the mass of Kepler- 
36b, and the impact velocity Vimp = Vesc, where Vesc = 

2G(Mj +M\) /(Rt +/^i) is the two-body escape velocity. In 
the disruptive model 2, Mt = 10 i.e. the upper end of the 
mass spectrum that is commonly refer to super-Earths, and 
the impact velocity vimp = 3vesc comparable to the local Ke¬ 
plerian speed. Our model 2 represents the largest scale of 
collisions that could happen to a close-in super-Earth with an 
atmosphere. Impacts onto a less massive super-Earth at the 
same speed lead to a more disruptive outcome. 

In our models, the target is initially composed of an iron 
core, a silicate mantle, and an H/He gaseous envelope (7.5 
wt %), while the impactor has the same layered structure ex¬ 
cept without an atmosphere. The mass ratio of an iron core to 
a silicate mantle is assumed to be 1: 2 for all planets. We sum¬ 
marize the masses of each species in the targets and impactors 
in Table 1 . The internal structure of all the target planets are 
resolved by a numerical mesh with Rt/AJ > 100, where 
is the target’s initial radius and AJ is the width of the smallest 
grid cells. 

Previously, Liu et al. (2015) have applied an early version 
of this scheme to investigate embryo impacts on gas giants. In 
those simulations, we approximated the internal structure of 
gas giants with a composite polytropic equation of state (EOS; 
Liu et al. 2013). We assume that the multi-species fiuid obeys 
the Dalton’s law, i.e. the total pressure of a mixture of gas 
equals the sum of the contributions of individual components. 
The current study focuses on collisions between smaller plan¬ 
ets with compositions and internal structures similar to those 
of the Earth or Neptune. Eor ’terrestrial’ types of planetary 
materials, we have incorporated the Tillotson EOS to model 
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TABLE 1 

Mass and energy budgets eor the low and high speed impacts. 


Timp 

(km/s) 

Iron’’ 

Rock" 

H/He^’ 

Ml 

Iron Rock 

Iron 

M 

Rock 

H/He 

Ut 

Ek,T 

Ui 

Ek,i 

U 


E^ 

rp rn 

10.96 

1.344 

2.663 

0.297 

0.334 

0.663 

1.676 

3.315 

0.208 

4.746 

0.441 

0.516 

2.092 

9.176 

9.442 

0.186 

0.215 

53.06 

3.115 

6.129 

0.746 

0.334 

0.662 

3.449 

6.289 

0.151 

21.890 

4.529 

0.407 

72.338 

31.039 

32.810 

0.240 

0.311 


^ Integral quantities M, U and Ek correspond to mass, internal energy and kinetic energy, respectively. Subscript T and I denote the target and the impactor. Integral quantities without 
additional subscripts are measured at the end point of our simulations. Mass data are in units of Earth mass and energy data are in units of 10^^ erg. 

c, d ]y[^gg Qf species iron, rock and H/He, respectively. 

^ The integral quantity is measured within the original radius of the target (cf. Figure 2 and 7). 

^ The integral quantity is measured within the Hill radius (2.62 x 10^® and 3.21 x 10^® cm, respectively). 



Fig. 1. — Snapshots of the giant impact simulation between a 4.3 super-Earth and an Earth-mass planet at the escape velocity Vesc- The density distribution 
across slices in the planets’ orbital plane with density range [10“^ , 2.5 x 10^] g cm“^ (represented by colors ranging from purple to white) before the impact 
(panel a), immediately after (b), at 1.56 hours (c), and at 18 hours after the impact (d). Panel e shows an enlarged view of the panel d. With a radius 2.62 x 10^^ cm, 
the dashed circle overplotted on the planet indicates its Hill sphere at a distance of 0.1 AU from a solar-mass star. The parent star is in the direction of the white 
arrow. 


iron cores and silicate mantles (see Appendix I of Melosh 
1989). We additionally model an H/He atmosphere (70 wt. % 
H 2 and 30 wt. % He) assuming a polytropic EOS {P (x p^) 
with an adiabatic index of 7 = 5/3. This simplified assump¬ 
tion can be validated for temperatures below the excitation 
temperature of the H 2 rotation band, 170K, and above H 2 dis¬ 
sociation temperature (^ 2000 K) but below that for the H and 
He atoms to be fully ionized. For the intermediate regime, 
our treatment underestimates the compressibility and heat ca¬ 
pacity of hydrogen, which may exaggerate an impact-driven 
atmospheric escape. However, the atmosphere of the target is 
generally not sufficiently massive to stall the impactor and to 
prevent the dissipation of its kinetic energy near the target’s 
core. We find that the computed turbulent mixing of material 
and the mass loss rate of the atmosphere are relatively insen¬ 
sitive to the specific H/He EOS model. 

Giant impacts are violent and impulsive events. A new 
quasi-hydrostatic equilibrium within the planet’s Hill’s radius 
is generally established within a dynamical timescale (within 


a day). During the transitory phase, we can neglect the effect 
of stellar heating on the inflated atmosphere (see §4). After 
the impact, the efficiency of insolation due to stellar irradia¬ 
tion on the planet’s surface and the efficiency of heat transfer 
in its interior determine the mass-loss rate from the inflated 
planet (Owen & Wu 2015). Here we must also consider that 
the giant impact has left behind a massive debris torus that 
can remain opaque for some time. 

Part of the analysis and visualization presented in Section 3 
is generated using the YT package (Turk et al. 2011) and the 
Visit software (Childs et al. 2012). 

3. RESULTS 
3.1. Low-speed model 1 

The compositional structure of a super-Earth after a low- 
speed collision with an Earth-mass impactor is shown in Fig¬ 
ure 1. The initial mass of the target is 4.3 M^. The im¬ 
pactor’s approaching speed is the escape velocity of the target- 
impactor system. Panels (a) - (c) show snapshots of density 
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Radius (cm) 

Fig. 2.— Each shaded area illustrates the mass fraction of each species as 
a function of radius at 18 hours after a low-speed giant impact. The grey, 
light coral, and blue regions represent iron, silicate rock and H/He gas, re¬ 
spectively. Three vertical dot-dashed lines indicate the initial compositional 
boundaries inside the target that has the three-layered structure. The dashed 
lines from top to bottom show the cumulative mass (yellow), the enclosed 
mass of silicate rock (red), iron (black), and H/He (blue) as a function of 
radius in units of Earth mass. 



Fig. 3.— Cutaway view showing a snapshot after the contact between the 
rocky layers of the two planets and before the merger of their iron cores. The 
blue, green and red colors represent partial densities of H/He, rock and iron 
species, whichever is the dominant species of a grid cell. For the most part 
of the target, the total density of a grid cell is close to the partial density of 
the dominant species because mixing is not severe. However, a small fraction 
of H/He (less than 5 wt. %) being mixed with rock is responsible for the 
formation of a low-density arced structure between compressed rocky layers. 
The cubic box has a width of 4 x 10^ cm. 


contours before, during, and at 1.56 hr after the impact. The 
snapshots in panels (d) and (e) are taken at 18 hours after the 
impact. At this time, the central region of the post-impact tar¬ 
get is mostly dynamically relaxed (see panel (d)), while the 
hot atmosphere extends well beyond the Hill radius and con¬ 
tinues to lose its mass via Roche lobe overflow (see also §4). 

Since the total orbital energy of the target-impactor system 
is nearly zero, this impact has an accretion efficiency close 


to but not quite one. The target gains mass after the impact 
(Table 1) but not all of it. The excess kinetic energy injected 
by the projectile drives an outflow as the planet approaches a 
hydrostatic equilibrium. As a result, the impact is not a per¬ 
fect merger. The total mass loss (i.e. the mass that is failed 
to be accreted by the target) at 18 hours after the impact is 
^ O.1M0. This value is expected to increase slightly due 
to the slowly decaying oscillations in the tenuous outer part 
of the atmosphere over a prolonged period of time. Accre¬ 
tion efficiency in low-velocity collisions can be substantially 
lower than it is for the head-on cases considered here, because 
off-axis giant impacts can attempt to accrete too much angu¬ 
lar momentum for a single planet or planet-satellite system to 
sustain. The so-called hit-and-run regime is relevant at veloc¬ 
ities intermediate to those considered here (Asphaug 2010). 

The gravitationally unbound mass, ^ O.IM 0 , is mostly 
composed of H/He gas plus a minute contribution from iron 
and silicate material (see Table 1). The impact leads to the 
immediate ejection of nearly one-third of the target’s initial 
H/He atmosphere. Most of the impactor’s mass is added to 
that of the target. Consequently, the volatile fraction of the 
target is significantly reduced. 

The enclosed mass as a function of radius inside the Hill 
sphere is plotted with a yellow dashed line in Figure 2. Most 
of the planetary mass is confined within its original radius, 
while the tenuous atmosphere spreads out by an order of mag¬ 
nitude, filling up the entire Hill sphere (see panel (e) of Figure 
1). The grey, pink, and blue regions in Figure 2 illustrate iron, 
rock and H/He mass fractions as functions of planet’s radius. 
The masses of the three components enclosed within a given 
spherical radius are overplotted with red (rock), black (iron), 
and blue (H/He) dashed lines, respectively. 

When the intruding planet enters the target’s envelope, a 
velocity shear is present at the interface between two species 
which triggers Kelvin-Helmholtz (K-H) instabilities. On the 
other hand, during the propagation of the impact-induced 
shock wave through the target, a denser fiuid is accelerated 
by a less dense one at their interface and Rayleigh-Taylor (R- 
T) instabilities start to develop. However, since the impactor 
merges with the target in less than half an hour, only small- 
scale instabilities arise during such a short time span (see, e.g. 
Agertz et al. 2007), which can only be resolved with an ex¬ 
tremely high resolution. Figure 3 shows a cutaway view of 
a snapshot when the impactor is about to merge with the tar¬ 
get. We plot the partial density, i.e. the density times the mass 
fraction of the dominant species of a grid cell, to illustrate 
the interior structure. We note that only a little mixing occurs 
when the target’s gaseous envelope gets crushed to its rocky 
layer, reducing the density of rocky material in a narrow re¬ 
gion between the target and the impactor (see Figure 3). As 
during the short impact phase, both K-H and R-T instabilities 
are unable to cause large scale mixing, we identify that the 
impact-driven turbulence is responsible for the global mixing 
observed during the late stages (see Figure 4). Because a large 
amount of kinetic energy is released upon the coalescence of 
the target and the impactor, fiuid motions become turbulent 
and a complex mixture of K-H and R-T instabilities can be 
triggered. The turbulent mixing, as well as hydrodynamic in¬ 
stabilities, destroy the original layered interior structure of the 
target - causing radial mixing between iron and rock near the 
center, and inducing a fraction of rock (and a much smaller 
fraction of iron) to diffuse into the H/He atmosphere (see Fig¬ 
ure 4). Nevertheless, figure 2 shows that more than two-thirds 
of the mass at the center still consists of iron species. Besides, 
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Fig. 4.— Mass fraction slices of each species through the orbital plane at 18 hr after the low-speed impact. A significant degree of mixing is established due to 
the turbulent mixing as well as hydrodynamic instabilities during the post-impact expansion phase. Dashed circles represent the the Hill sphere with a radius of 
2.62 X 10^^ cm. 



Fig. 5.— Compositional gradients inside the target planets after head-on 
collisions. The low-speed and high-speed impacts are represented by solid 
and dashed lines, respectively. Mass fractions of iron (black), rocky material 
(red), and heavy elements (blue) as functions of the radius are shown with 
labels on the right. The compositional gradient Vz (green) increases with 
radius in the low-speed case. The small value of Vz at radii indicates that the 
elements are well-mixed and homogenized. Fluctuation in Vz reflects that 
the planet has not yet reached a relaxed state. 

the iron mass fraction falls off rapidly with the radius, indi¬ 
cating that the planetary core survives from the impact and its 
mass grows in a coalescent manner. 

In contrast, the rocky material is dredged up to the top of the 
atmosphere due to the turbulent motions driven by the impact. 
As shown in Figure 5, we note that the target’s interior, i.e. 
within 6x10^ cm, has a steep and positive compositional gra¬ 
dient, Vz = dlogZ/dlogF, which serves to stabilize thermal 
convection in the mantle. At least right after the giant impacts, 
a head-on collision develops a hot and inhomogeneous inte¬ 
rior. Heat transport deep inside a post-impact planet may be 
suppressed by the effect of double diffusion convection (e.g. 
Rosenblum et al. 2011; Leconte & Chabrier 2012). Such an 
inefficient heat transfer mechanism would also prolong the 
presence of a magma ocean and delay the differentiation of 
iron material from a silicate-rich mantle. Off-axis collisions 
would lead to more efficient mixing, but also a less efficient 
kinetic energy contribution. 

We summarize the energy budget in Table 1 . Note that the 


planet’s internal energy is almost doubled after the impact and 
it is larger than the sum of the internal and kinetic energies of 
the target and the impactor prior to the impact. This differ¬ 
ence arises because the collision is accretionary, i.e. a large 
amount of the impactor’s mass is delivered to the interior of 
the target. As a result, not only most of the kinetic energy but 
also some portion of the gravitational potential energy is con¬ 
verted into the post-impact internal energy. Therefore, one 
would expect that the temperature of the planet’s interior to 
increase substantially after the impact. 

3.2. High-speed model 2 

The consequences of a high-speed impact is similar to that 
of the low-speed counterpart in terms of shock propagation 
throughout the target after the impact, mass loss via Roche 
lobe overflow, and gradually decaying oscillations inside the 
planet (Figure 6). Despite those general similarities, the high¬ 
speed impact differs from the low-speed one in three respects. 

First, the high-speed impact is erosive. The mass of the 
target planet at 21.5 hr after the impact is about 9.9M0, i.e., 
O.IM 0 is stripped off from the target by the impactor. Thus, 
the total mass loss is ^ I.IM 0 . Besides, a significant amount 
of mass is levitated at the top of the atmosphere (see Figure 
7). A breakdown of the mass budget in the Hill sphere (i.e., 
3.21 X10^^ cm in this case) shows that the impactor success¬ 
fully deposits all of its iron material into the target, while most 
of the impactor’s rocky material is ejected (see also Table 1). 
And the target retains about 1 /5 of its original atmosphere 
after such an energetic impact. 

Second, silicate rock material become the dominant species 
in the planet’s atmosphere (Figure 7 and 8). As a result, the 
compositional gradient is nearly zero out to the Hill’s radius, 
indicating a homogeneous distribution of the heavy material 
inside the bulk of the target (see Figure 5). In contrast to 
the low-velocity model, model 1, the negligible compositional 
gradient cannot suppress convection in the interior of the tar¬ 
get planet. Nevertheless, iron material remains concentrated 
near the center of the planet as in the case of the low-speed 
model, model 1. 

Third, the planet becomes less bound gravitationally. The 
mass contained within the target’s original size is much 
smaller (see Figure 7). In terms of the energy budget, most 
of the kinetic energy is carried away by the unbound mass 
and only a fraction of it is directly converted to the planet’s 
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Fig. 6. — Snapshots of the giant impact simulation between a 10 super-Earth and an Earth-mass planet at 3 Vesc- The colormaps and symbols are same as 
in Figure 1. Panels (a)-(d) are snapshots taken at the start of the impact, at 15 minutes, at 1.5 hr and at 21.5 hr after the impact, respectively. Panel (e) shows an 
enlarged view of the panel (d). The Hill sphere illustrated by the dashed circle has a radius of 3.21 x 10^^ cm. 

neous outcomes depending on energy of collision, one at low 
velocity and the other at approximately the local Keplerian ve¬ 
locity, three times the planets’ escape velocity. Whereas the 
initial gaseous envelope of a target super-Earth is mostly re¬ 
tained during the low-velocity impact, it is severely depleted 
shortly after the high-velocity impact. In this context, we sug¬ 
gest that giant impacts of varying energy can effectively de¬ 
volatilize super-Earths/mini-Neptunes, and can diversify the 
compositions and interior structures of these mid-sized ex¬ 
trasolar planets, perhaps as occurred in our terrestrial system 
(see, e.g. Asphaug 2014). 

The interior and atmosphere of a post-impact planet can 
achieve a thermal equilibrium state after a long-term evolu¬ 
tion, which is difficult and not cost-effective to study by per¬ 
forming hydrodynamic simulations. Alternatively, the radia¬ 
tive cooling timescale (rj-ad) can be estimated from the rate of 
temperature change in an atmospheric layer via the outgoing 
infrared radiation: Trad ^ PC^/{Ag(7T^), where P and T are 
the pressure and temperature near the photosphere obtained 
from the end state of hydrodynamic simulations, g is the sur¬ 
face gravity, and Cp is the specific heat capacity at a constant 
pressure, and a is the Stefan-Boltzmann constant. Under the 
intense radiation from a 1 Mq central star, an isothermal layer 
develops in the planet’s residual upper atmosphere. 

At the current location of Kepler 36b and c, the equilibrium 
temperature due to stellar irradiation is ^ 900 K. Adopting 
EOS of an ideal gas for H/He and the gas opacity of Ereedman 
et al. (2008), we find that the H/He gas within the Hill sphere 
is optically thick. At the photosphere, we used P ^ 2gl?>n 
and n is the opacity to estimate Trad to be larger than several 
days. This estimate justifies, a posteriori the neglect of stellar 
irradiation in our impact simulations (see §2). Since this cool¬ 
ing time scale is much longer than the dynamical time scale, 
we can also assume that after a giant impact, a heated planet 


J ^ iron ; KOCJC ; LL/Me ; _-—I 10 
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Fig. 7.— Mass fraction area plot and enclosed mass line plot at 21.5 hours 
after the high-speed impact. The symbols have the same meanings as in Fig¬ 
ure 2. 

internal energy. However, the gravitational potential energy 
released by the sedimentation of rock material in the atmo¬ 
sphere becomes an extra heat source for the target interior. 
Assuming that silicates in the H/He atmosphere establish a 
local thermal equilibrium with the turbulent gas, they con¬ 
dense into grains in cooler outer regions with r > 10^^ cm. 
These silicate condensates grow through collisional coagula¬ 
tion (e.g. Podolak 2003). Relatively large silicate grains are 
expected to settle toward the planetary surface and sublimate 
along the way. 

4. DISCUSSION AND CONCLUSIONS 

The results presented in the previous sections show the 
planet’s structure immediately after head-on giant impacts. 
Models 1 and 2 illustrate the possibility of diverse instanta¬ 
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Iron Silicate H/He 



Fig. 8. — The mass fraction slices of each species through the orbital plane at 21.5 hours after the high-speed impact. Mixing becomes less efficient because 
most of the H/He envelope has been blown away by the impactor. Dashed circles represent the Hill sphere with a radius of 3.21 x 10^^ cm. 


evolves adiabatically into a hydrostatic equilibrium. 

During the post-impact phase, the planet’s atmosphere and 
interior contract as the planet cools. Above the photosphere, a 
highly-irradiated planet has an isothermal layer. We consider 
that an inflated planet after an impact radiates away heat from 
the photosphere. Owen & Wu (2015) estimated the Kelvin- 
Helmholtz timescale of a planet (3, 5, and IOM 0 ) with the 
equilibrium temperature of 900 K as a function of planetary 
radius in unit of the Bondi radius and envelope mass fraction 
based on a stellar evolution code, i.e., MESA (Paxton et al. 
2011). We applied their results of 5 and IOM 0 shown in Fig¬ 
ure 2 in Owen & Wu (2015) to our targets after giant impacts 
(^ 4.99IM 0 and 9.738M0 with the equilibrium temperature 
of-900 K). 

Following results of Owen & Wu (2015), we estimate the 
timescale of the Kelvin-Helmholtz contraction to be — 1 Myr 
for the low-speed case and <10 kyr for the high-speed one. 
A protracted state of a hot and inflated atmosphere enhances 
its mass loss via a Parker wind. Based on Eq.(16) in Owen & 
Wu (2015), we estimate that the amount of mass loss for post¬ 
impact planets via a Parker wind for the two cases would be 

— 88% and — 80% of its envelope mass right after the impact. 

As the planet cools down, silicate material sediments from 

its atmosphere and contracts within its interior. The release of 
gravitational energy provides a source of internal heat which 
prolongs the inflationary state of the planet’s atmosphere and 
enhances mass loss via a Parker wind. The extended planetary 
atmosphere also intercepts a greater fraction of the stellar X- 
ray and UV irradiation and increases the rate of mass loss 
through photoevaporation. 

We estimate the mass loss rate of an extended atmosphere 
via stellar XUV irradiation after giant impacts. Adopting an¬ 
alytical descriptions of a XUV energy flux from solar-type 
stars (Ribas et al. 2005), the incident XUV luminosity at the 
Hill radius is 1.0 x 10“^ Lq and 1.5 x 10“^ L 0 in the low-speed 
and high-speed model, respectively. Given that the conversion 
efficiency from XUV photons to the kinetic energy of bulk at¬ 
mospheric outflow is 10% as a flducial value (Yelle 2004), the 
mass loss rate from the Roche lobe via stellar XUV irradia¬ 
tion right after giant impacts is estimated to be — 3M0 Myr“^ 
for the former and — 2 M 0 Myr“^ for the latter. If we consider 
a typical decay timescale of a XUV flux for Sun-like stars is 

— 0. IGyr, the planet in the low-speed model may lose the en¬ 
tire envelope because its atmosphere shrinks in a few Myr. On 


the other hands, the planet in the high-speed model can con¬ 
tract more quickly (< 10kyr), but it is unlikely for its tenuous 
post-impact atmosphere to survive a subsequent mass loss via 
stellar irradiation, either. Note that the stellar tidal held is cru¬ 
cial in this context because it can continuously remove plan¬ 
etary outer atmosphere that is beyond the Hill sphere. For 
giant impacts happen further away from the host star, tidal 
stripping becomes inefficient. And the outcome of a giant im¬ 
pact is determined by the mass ratio, impact speed and impact 
angle (Asphaug 2010; Feinhardt & Stewart 2012; Stewart & 
Feinhardt 2012). 

To summarize, severe giant impacts can signiflcantly de¬ 
volatilize close-in super-Earths and mini-Neptunes. The com¬ 
positional dichotomy between Kepler-36 b and c can be ex¬ 
plained by their distinct impact histories along with the for¬ 
mation of the closely packed system, i.e. Kepler-36b experi¬ 
enced substantial giant impacts, and in the meanwhile Kepler- 
36c survived from being heavily bombarded. In addition, we 
speculate that giant impacts may have been imprinted in the 
large dispersion in the mass-radius relationship of close-in 
sub-Neptune-sized planets, as giant impacts occur stochasti¬ 
cally and can diversify planetary interior and atmosphere oth¬ 
erwise well constrained. 

In our solar system, there may be evidence for a similar 
though more subtle dichotomy. Despite their similar masses, 
radii, and compositions, thermal evolution models of Uranus 
and Neptune suggest that right after their formation, Nep¬ 
tune may have been relatively luminous but Uranus relatively 
faint (e.g. Hubbard & Macfarlane 1980; Fortney et al. 2011). 
Stevenson (1986) suggested that this dichotomy may be ac¬ 
counted for if a violent head-on collision yielded a hot and 
homogeneous interior of Neptune, whereas an oblique col¬ 
lision caused a tilted Uranus with a stably-stratifled interior. 
The results presented here imply that the former scenario for 
Neptune may be possible. 
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